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ABSTRACT 

We present a model for chemical evolution of damped Lya systems considering production of "metals" 
by SNe II and infall associated with hierarchical structure formation. The growth of metallicity in these 
systems is a reflection of the competition between astration and infall. The apparent late turn-on of 
these systems is due to the late cut-off of infall. The wide range in [Fe/H] at a given redshift is explained 
by the range of the times for onset of star formation and the range of the times for infall cessation in 
different systems. The observed lower bound of [Fe/H] « —3 follows from the very rapid initial rise 
of [Fe/H] subsequent to onset of star formation. To reach [Fe/H] « —3 from a metal-free initial state 
requires only 30 Myr so that the probability of observing lower [Fe/H] values is very small. 

Subject headings: galaxies: abundances — galaxies: evolution — intergalactic medium 



L INTRODUCTION 

We discuss chemical evolution of damped Lya (DLA) 
systems based on hierarchical structure formation. Ob- 
servations (Lu et al. 1996; Prochaska & Wolfe 2000, 
2002; Prochaska et al. 2003) show there is a baseline 
enrichment of [Fe/H] = log(Fe/H) - log(Fe/H)o « -3 
for DLA systems. Of 96 systems over the redshift range 
0.5 < z < 5, the lowest observed [Fe/H] is —3.13 although 
the newer observations could have detected much lower 
values (Prochaska et al. 2003). Furthermore, there is a 
large dispersion in [Fe/H] at any fixed z. 

In earlier discussion (Wasserburg & Qian 2000b), we 
treated DLA systems as closed systems formed from the 
intergalactic medium (IGM) with a fixed initial inven- 
tory of "metals" (the prompt inventory) corresponding to 
[Fe/H]p « —3 and assumed that subsequent Fe enrich- 
ment of each system was governed by a constant Fe pro- 
duction rate per H atom in the gas of PFe/(H). Element 
production started when a system was formed at time t* 
after big bang. At time t > t* , the number ratio (Fe/H) 
of Fe atoms to H atoms in the gas evolves according to 



d(Fe/H) 
di 



Pfc 
(H)- 



(1) 



We only treat Fe contributions from Type II supernovae 
(SNe II) in all our considerations as SNe la turn on at later 
times. Integration of equation (1) gives 



Pi 



(Fe/H) = (Fe/H)p + ^(i-t*) 



(2) 



As SNe II contributed ^ (Fe/H)Q/3 over the period of 
(t — t*) ^ 10 Gyr prior to solar system formation, Ppc/ (H) 
is estimated as ~ (Fe/H)g/(30 Gyr). 

In this model, the baseline Fe enrichment is explained 
by the prompt inventory and the dispersion in [Fe/H] at a 
fixed z by the range of t* . It was found that most systems 
at any given z have t* close to the age t{z) of the universe 



at this z. This implies that the "turn-on" of DLA systems 
occurs rather long after the big bang. Further, essentially 
all the data lie below the upper bound for [Fe/H] corre- 
sponding to t* — 0. The above model appears to provide 
a reasonable description of the data on [Fe/H] for DLA 
systems. This is rather remarkable considering that the 
model ignores infall. According to hierarchical structure 
formation, infall is essential to formation of DLA systems. 
Here we examine chemical evolution of these systems by 
including infall. We wish to gain some insights into what 
causes most DLA systems at a given z to have t* close to 
t{z) in the closed-system model. We also address whether 
the prompt inventory is needed to explain the baseline Fe 
enrichment for DLA systems. 

2. CHEMICAL EVOLUTION WITH INFALL 

Consider a system of gas and stars with infall of pri- 
mordial metal-free gas. The equations for evolution of the 
numbers of Fe and H atoms [(Fe) and (H)] in the gas are 



d(Fe) 
dt 

m 

dt 



Fc 



(fc/h; 



rf(H)a 



dt 



d(H)as , d(H)i, 



dt 



dt 



(3) 
(4) 



where Ppc is the Fe production rate of SNe II, d(H)s,s/dt < 
is the astration rate, and d{H)ia/dt > is the infall rate. 
When ]d(H)as/rfi| is small compared to d(H)in/dt, (H) is 
governed by infall. Equations (3) and (4) give (Qian & 
Wasserburg 2003) 



d(Fc/H) 
dt 



(H) 



1 d(H)i; 

(H) dt 



■(Fe/H). 



(5) 



Formal integration of equation (5) gives 



(Fe/H)^^t- 



1 d(H)i, 



■{Fc/R)dt'. 



(H) dt' 
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With the integral term interpreted as Pfo^*/(H), this re- 
sembles equation (2) for (Fe/H)p = 0. So the effect of a 
sustained large infall rate is similar to late "turn-on." 

In general, equations (4) and (5) must be solved to- 
gether. For simplicity, we first consider the case where 
d{B.)in/dt > \d{R)^/dt\ so (H) « (H);^. We then have 



dZp^ 
dt 



(7) 



where Zp, = (Fe/H)/(Fe/H)0, Ape = PFe/[(Fe/H)Q(H)], 
and Ain = d \ii(K) in /dt. If dZ^c/dt is negligible, a quasi- 
steady state is achieved and .^Fe(0 approximately assumes 
the value Z^^^{t) = XFe{t)/Xin{t). With onset of Fe pro- 
duction at to, we have 



ZFeit) 



zQf (f) 



Z^/S(io)exp 



Xin{t')dt' 



to 



(8) 



the 



for Ai„(t) > \d\nZ§^^/dt\. At t > to + [Ain (to)]" 
exponential term is negligible and Zpeit) w Zp^^{t). 

We consider a baryonic system that is formed through 
infall into the potential well of a dark matter halo. The 
mass M of the halo grows according to hierarchical struc- 
ture formation. We assume that baryonic and dark matter 
are fed into the halo at a fixed mass ratio. The infall rate 
is then 

, dhiM dlnMdz 

Ai„(t) = — — = -j— -. (9) 



dt 



dz dt 



With t 



17(1-1- z) Gyr at z > 0.5, equation (9) gives 



Ain(t) 



4£ 




5/3 


dlnM 


Gyr 




) 


dz 



(10) 



From Figure 6 in Barkana & Locb (2001), we find that 
\dlnM/dz\ ~ 2.3 (within a factor of 2) for 0.5 < z < 5. 

As an example, we take Xinit) = (0.1 Gyr)~^(Gyr/t)^/^ 
and Afo = (30 Gyr)-i for t > (i.e., to = 0). We nu- 
merically integrate equation (7) and show the evolution of 
[Fc/H] = logZpc in Figure 1. The quasi-steady state value 

Z^^^{t) = AFo/Ai„(t) = (i/Gyr)5/3/300 is a good approx- 
imation to the exact solution. For the case of evolution 
from a metal- free initial state without infall, Zp^it) = Xpet 
(the case of t* = in §1). In comparison with this, Fe en- 
richment in the case of infall is suppressed by a factor 
of « Ain(i)i « 10(Gyr/t)^/''*, which ranges from 4.5 for 
t = 3.3 Gyr (z = 2) to 8.9 for t = 1.2 Gyr (z = 5). This 
would explain the late "turn-on" requiring t* to be close 
to t{z) in the closed-system model. The data on [Fe/H] 
for 96 DLA systems (Prochaska ct al. 2003) arc shown in 
Figure 1. It can be seen that the solid curve for the case of 
infall passes through the median of the body of the data. 
This curve has a slope of d[Fe/B.]/dz « -1.1/(1 -|- z). For 
2^2 where most of the data lie, this slope is in good 
agreement with the estimate of d[Fc/}i]/dz « —0.26 given 
by Prochaska et al. (2003) for the growth of the mean 
[Fe/H] as a function of z. However, the above result does 
not explain the wide range in [Fe/H] at a given z. 

3. THE DISPERSION IN [FE/H] 

We now turn to the wide range in [Fe/H] for DLA sys- 
tems at a given z. We assume a constant Fe production 



rate Apo per H atom in the gas for t > t{). The infall rate 
Ain is estimated based on hierarchical structure formation 
and we now allow the possibility that infall may cease at 
a time tiQ. We assume that the infall rate greatly exceeds 
the astration rate at i < tie. Under these assumptions, to 
and tic determine [Fe/H] in a baryonic system at time t 
(see eqs. [8] and [11]). 

First consider the case where to = and infall ceases at 
tic- The evolution of Zpe is the same as for the case of 
infall discussed in §2 imtil t = tic. At t > tic, equation 
(7) reduces to dZpe/dt = Ape, which gives 



Zveit) = ■^Fe(tlc) + Ape(t — tic)- 



(11) 



Solutions for tic = 0, 1.5 Gyr (zic = 4), and 4.3 Gyr 
{zic = 1-5) are shown in Figure 2. In all cases of tjc > 0, 
the evolution of [Fe/H] for t < tic is along the solid curve 
labeled "continuous infall." As illustrated by points A, 
B, and C at z = 3 in Figure 2, a wide range of [Fe/H] 
bounded by the cases of no infall (tjc = 0) and continuous 
infall can be produced at a given z for different tic values. 

Next consider the case where infall is continuous but 
to may vary. Subsequent to onset of Fe production at to, 
Zpe grows rapidly and approaches Z^^^ on a timescale of 
~ Ain(to)~"'^ (see eq. [8]). This is shown in Figure 1 for 
to = 1.5 Gyr {zq = 4), 2.1 Gyr {zq = 3), and 3.3 Gyr 
{zo = 2). For example, [Fe/H] in baryonic systems formed 
at to = 2.1 Gyr first evolves rapidly, crossing [Fe/H]=-3, 
and then grows more slowly toward the quasi-steay state 
solution. Thus, values of [Fc/H] below those for the case 
of to = may be populated with baryonic systems that 
have different to values. Then the lower range in [Fe/H] 
at a given z can also be explained by the infall model. 

The observed lower bound of [Fe/H] w —3 for DLA sys- 
tems requires attention. With infall, the growth of Zpc 
is determined by competition between Ape and Ain-^Pe- 
As Zpe(to) = 0, Ape governs the initial growth of Zpe 
so long as Apo/Ajn ^ Zpc- We note that Apc/Ain = 
(t/Gyr)5/3/300 > lO'^ for t > 0.49 Gyr (z < 9.6). So for 
systems with to > 1-2 Gyr (zo < 5), the growth of Zpe up 
to 10~^ is essentially governed by Ape = (30 Gyr)~^. These 
systems achieve Zpc = 10~^ in ~ 30 Myr and then ap- 
proach Zpc « Zp^^ > 10^^ on a timescale of ^ Ain(to)^^ 
(see eq. [8] and Fig. 1). Thus, all baryonic systems 
will reach Zp^ = 10~^ or greater in a very short time at 
z < 5 and the probability for observing DLA systems with 
[Fe/H] < —3 is quite small. 

The ranges of to and tjc in the infall model to explain 
the dispersion in [Fe/H] at a given z are analogous to the 
range of t* in the closed-system model. However, uniform 
growth of [Fe/H] occurs only at t > tic in the infall model. 
It is the dilution of gas by fresh infalling baryonic mate- 
rial that retards [Fe/H] from uniform growth although Fe 
production and astration has been going on since to- To 
relate to and tic for a baryonic system to properties of 
the halo hosting the system, we follow the discussion in 
Barkana & Loeb (2001). The condition for a halo of mass 
M associated with an na density fluctuation (an na halo) 
to collapse at redshift z is 



1.33(1 -I- z) « na{M), 



(12) 



3 



where cr(M) is shown in Figure 5 of Barkana & Loeb 
(2001). Equation (12) gives the z values at which a fixed 
mass M is reached by halos with different n values. 

A IfT halo today has M w 10^^ Mq and should contain 
« IQ-'^^ Mq of baryonic matter. This mass is far greater 
than found for a typical galaxy. Thus, we may consider 
that infall into individual protogalaxies ceases when the 
accreted baryonic matter reaches ^ 10^°-10^^ Mq in a 
halo of ~ 10^^-10^^ Mq. In the hierarchical model, a halo 
of larger masses (e.g., ^ IQ-'^^ M0) may consist of a num- 
ber of smaller collapsed halos instead of a single object. 
This also applies to the associated baryonic systems. We 
assume that infall ceases at tic when a halo reaches a 
mass of Mic = 10" [cr(Mic) = 3.5]. This occurs at 
redshifts z\c = 4 and 1.5 for 1.9(7 and 0.95(t halos, respec- 
tively. These redshifts correspond to the t\c in Figure 2 
and seem reasonable. For > 4(7 halos, Mjc = 10^^ Mq is 
reached at 2:1c > 9.5 (tic < 0.5 Gyr). The evolution of 
[Fe/H] at < 5 (t > 1.2 Gyr) for baryonic systems inside 
those high a halos is close to the case of tic = 0. Similar 
results are found for Mic ~ lO^-lO^^ 

Now consider to, which may be taken as the onset of star 
formation. Following the discussion in Barkana & Loeb 
(2001), the first stars appear to have formed in ^ 3(7 to 4(7 
halos of mass M ~ 10^ Mq [(7(M) k. 10] at z - 20-30. For 
these halos, to is only ~ 0.1 Gyr. Formation of stars at 
2 < 20 appears to require a minimum halo mass of ~ 10^- 
10^ Mq. Wc assume that star formation starts at to when 
a halo reaches a mass of Mq = 5 x 10^ [(7(Mo) = 6]. 
This occurs at redshifts zq = 4, 3, and 2 for l.lcr, 0.89(7, 
and 0.67(7 halos, respectively. These redshifts correspond 
to the to in Figure 1 and again seem reasonable. Similar 
results are found for Mq lO^-lO^Af©. 

It follows that at a given z, baryonic systems inside dif- 
ferent halos are in different stages of evolution. Consider 
an example using Mo = 5 x 10^ Mq and Mic = lO" Mq. 
At 2 = 2.6, star formation just starts in a 0.8(7 halo 
([Fc/H] ^ ^3), Fc production competes with infall in 
a 1.1a halo ([Fe/H] = —1.9, see Fig. 1), and uniform 
growth of [Fe/H] has been going on since infall cessation 
at zic = 4 in a 1.9(7 halo ([Fc/H] = -1.4, see Fig. 2). For 
a sample of DLA systems at z — 2.6, systems with —3 < 
[Fe/H] < -1.9, -1.9 < [Fe/H] < -1.4, and [Fe/H] > -1.4 
are then associated with na halos with 0.8 < n < 1.1, 
1.1 < n < 1.9, and n > 1.9, respectively. Statistically 
speaking, the fraction of halos more evolved than an na 
halo is F{n) = y'2/7r J°° exp(— x^/2)(ia;. Thus, we expect 
the occurrences of -3 < [Fe/H] < -1.9, -1.9 < [Fe/H] < 
— 1.4, and [Fc/H] > —1.4 at z = 2.6 to be in the ratios of 
[F(0.8)-i^(l.l)] : [i^(l.l)-F(1.9)] : F(1.9) « 0.7 : 1 : 0.3. 
Of 19 DLA systems at 2.4 < z < 2.7, the numbers of 
systems in the above [Fe/H] intervals are 5, 10, and 4, 
consistent with expectation from our model. 

4. DISCUSSION AND CONCLUSIONS 

We have treated chemical evolution of DLA systems con- 
sidering astration and infall. It is assumed that the H in 
the gas is initially controlled by infall of primordial metal- 
free baryonic matter and the Fe production rate per H 
atom in the gas is constant subsequent to onset of star for- 
mation in a system. With the infall rate estimated from 
the standard scenario of hierarchical structure formation, 



this model yields an explanation for the data on [Fc/H] in 
DLA systems. It is shown that the slow growth of [Fo/H] 
with decreasing z is a direct consequence of competition 
between enrichment by Fe production and dilution by in- 
fall. It is argued that the upper range of [Fe/H] at a given 
z results from the different times at which individual halos 
associated with protogalaxies reach a mass of ~ 10^^ Mq 
and infall ceases or greatly diminishes. It is also argued 
that the lower range of [Fc/H] results from the different 
times at which individual halos reach a minimum mass of 
~ 10*-10^ Mq required to initiate astration. Because the 
initial growth of [Fe/H] up to —3 is very rapid subsequent 
to onset of astration, the probability of observing DLA 
systems with [Fe/H] < —3 is very small. The approach 
presented here appears to be a reasonably quantitative 
description of chemical evolution of DLA systems that is 
compatible with the paradigm of hierarchical structure for- 
mation. Of the available data, only three points lie outside 
the upper bound of the model and these are still within a 
factor of 2 of this bomid. 

We assumed the astration rate is small compared to the 
infall rate during the infall phase. If we assume a con- 
stant astration rate a per H atom in the gas consistent 
with the Fe production rate, then equation (4) reduces to 
(i(H)/fit = —a(H) + d(H)in/dt, which can be solved explic- 
itly. As long as d(R)in/dt > a{H), (H) increases. When 
infall ceases, (H) reaches the maximum value (H)max and 
the star formation rate also reaches the maximum value 
'^(H)max- Subsequently, the gas will be depleted by star 
formation on a timescale of ~ a~^. When a protogalaxy 
would reach its maximum star formation rate depends on 
when its associated halo reaches a mass of ~ 10^^ Mq and 
infall ceases. This occurs at z 1.6 and 4.3 for la and 2a 
halos, respectively. The extent to which the approach out- 
lined here can be used to quantitatively explain the cosmic 
star formation history remains to be explored. 

In earlier works (Wasscrburg & Qian 2000a; Qian & 
Wasserburg 2002), we proposed that very massive (;> 
100 Mq) stars (VMSs) produced the elements from C to 
the Fe group to explain the observed jump in the abun- 
dances of heavy r-proccss elements (Ba and above) at 
[Fe/H] « —3. The VMSs produced no heavy r-elements 
but dominated chemical evolution at [Fe/H] < — 3. This 
evolution resulted in a prompt inventory of "metals" cor- 
responding to [Fe/H] « —3 in the IGM. Cessation of 
VMS activities and rapid occurrence of a hypothesized 
SN II source for heavy r-elements then led to an expla- 
nation for the observed jump in abundances of these el- 
ements at [Fe/H] « —3. The dominance of VMS ac- 
tivities at [Fe/H] < —3 and transition to normal astra- 
tion at [Fe/H] « —3 are in accord with the theoretical 
study of Bromm et al. (2001). They showed there is a 
critical metallicity corresponding to —4 < [Fe/H] < —3 
above which normal star formation should occur. Oh ct 
al. (2002) showed that the number of VMSs required to 
produce [Fe/H] « — 3 would also provide enough photons 
to reionize the IGM. Qian, Sargent, & Wasserburg (2002) 
showed that the abundances of C, O, and Si in Lya forests 
appear to be consistent with the corresponding prompt 
inventory. The observed lower bound of [Fe/H] fa —3 
for DLA systems was interpreted by Wasserburg & Qian 
(2000b) as strong evidence in support of the prompt in- 
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ventory model. Thus, there was a general accord between 
very different approaches. 

There are now two considerations that require a re- 
evaluation of the prompt inventory model. Qian &: Wasser- 
burg (2003) proposed that the possible sources for heavy 
r-elements arc accretion-induced collapse (AIC) of white 
dwarfs and low-mass (~ 8-10 Mq) SNe II, with the most 
likely source being AIC events. If AIC events instead of 
SNe II are the dominant source for heavy r-elements, the 
jump in abundances of these elements at [Fe/H] w —3 
can be attributed to the delay in turn-on of these events 
relative to SNe II without invoking VMSs. We showed 
here that the growth of [Fe/H] to —3 from a metal-free 
initial state is very rapid (in ~ 30 Myr). So the proba- 
bility of finding DLA systems with [Fe/H] < — 3 is very 
low and the observed lower bound of [Fe/H] —3 can 
no longer be considered as an argument in favor of the 
prompt inventory model. A prompt inventory correspond- 
ing to [Fe/H] « —3 is thus not required to explain the data 
on DLA systems but is also not in conflict with the inter- 



pretation of these data using the infall model. However, if 
there is no prompt inventory, a question remains regarding 
the approximately constant metallicity of the IGM over 
1.5 < ^ < 5.5 (Songaila 2001; Pettini et al. 2003). At 
present, it is not obvious how this can be explained by the 
infall model. Without VMSs associated with production 
of the prompt inventory, the source for reionization may 
also be a problem. 
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galax;ies of reasons. We thank Jason Prochaska for mak- 
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lus. This work was supported in part by DOE grants DE- 
FG02-87ER40328, DE-FG02-00ER41149 (Y. Z. Q.) and 
DE-FG03-88ER13851 (G. J. W.), Caltech Division Con- 
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Fig. 1. — Data (squares) on [Fe/H] for 96 DLA systems (Prochaska et al. 2003). The solid curve labeled "infall" shows the evolution of 
[Fe/H] for a system with infall and a constant Fe production rate Ape per H atom in the gas since the big bang. The thin dashed curve is the 
quasi-steady state approximation. The solid curve labeled "no infall" is for a system with a constant Xpg since the big bang but with no infall. 
The short-dashed curve is for a system in which infall starts at the big bang but star formation starts at to = 1.5 Gyr (zq = 4), the dot-dashed 
curve is for to = 2.1 Gyr (zo = 3), and the long-dashed curve is for to = 3.3 Gyr (zo = 2). Assuming a halo mass of Mo = 5 X 10** Mq at the 
onset of star formation, zo = 4, 3, and 2 correspond to l.lir, 0.9a, and 0.7o- density fluctuations, respectively. Due to very rapid initial rise of 
[Fe/H], the probability of sampling DLA systems with [Fe/H] < —3 is very small. The prompt inventory proposed earlier (horizontal dashed 
line) is shown for reference. 

Fig. 2. — Data and the two solid curves same as in Figure 1. The dot-dashed curve shows the evolution of [Fe/H] at t > tic for a system 
in which infall starts at the big bang but ceases at tic = 1-5 Gyr (zic = 4) and the dashed curve is for tic = 4.3 Gyr (zic = 1-5). The 
overall evolution of [Fe/H] for the system with tic = 1-5 Gyr is indicated by the arrows. In all cases, star formation starts at the big bang 
(to = 0). Points A, B, and C at z = 3 represent [Fe/H] values sampled by DLA systems with different tic values. Assuming a halo mass 
of Mic = lO^'^ Mq (Ri lO'^" Mq in baryonic matter) at the time of infall cessation, zic = 4 and 1.5 correspond to 1.9a and 0.95it density 
fluctuations, respectively. 



-0.5 



-1 



-1.5 



5 



-2 



-2.5 - 



-3 



-3.5 




0.5 



1.5 



2.5 



3.5 



4.5 



